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Abstract. We describe how star formation is expected to proceed in tlie early metal-free Universe, focusing on the very 
first generations of stars. We then discuss how the star formation process may change as the effects of metallicity, external 
radiative feedback, and magnetic and turbulent support of the gas become more important. The very first stars (Pop III. 1) have 
relatively simple initial conditions set by cosmology and the cooling properties of primordial gas. We describe the evolution 
of these stars as they grow in mass by accretion from their surrounding gas cores and how the accretion process is affected and 
eventually terminated by radiative feedback processes, especially H II region expansion and disk photoevaporation. The ability 
of the protostar and its disk to generate dynamically important magnetic fields is reviewed and their effects discussed. Pop 
III. I star formation is likely to produce massive (~ 100 — 2OOM0) stars that then influence their surroundings via ionization, 
stellar winds, and supernovae. These processes heat, ionize and metal-enrich the gas, thus altering the initial conditions for 
the next generation of star formation. Stars formed from gas that has been altered significantly by radiative and/or mechanical 
feedback, but not by metal enrichment (Pop III.2) are expected to have significantly smaller masses than Pop III. I stars because 
of more efficient cooling from enhanced HD production. Stars formed from gas that is metal-enriched to levels that affect the 
dynamics of the collapse (the first Pop II stars) are also expected to have relatively low masses. We briefly compare the above 
star formation scenarios to what is known about present-day star formation. 
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INTRODUCTION 

The first generation of stars is expected to have played 
a significant role in reionizing the Universe, for which 
observational constraints can be derived from CMB po- 
larization 1 1] and future high redshift 21 cm HI observa- 
tions [2]. The first stars polluted their pristine surround- 
ings with metals and may have contributed to the ob- 
served metal abundances of Galactic halo stars yl] and 
the Ly-a forest |4, 5]. Light from the first stars may con- 
tribute to the observed NIR background intensity, e.g. 
la, ED, and its fluctuations fl (however, see ^). If mas- 
sive, the deaths of the first stars may be observable as 
supernovae 1101 or gamma-ray bursts |J_1]. The first stars 
thus set the stage for galaxy formation, for which obser- 
vational constraints now exist at redshifts 8 to 10 [121], 
and supermassive black hole formation and growth, with 
'^billion solar mass black holes being seen out to redshift 
6.4 tllH. 

In discussing the formation of the first stars, the terms 
"first stars" and "Population III stars" are often used in- 
terchangeably, but this can lead to confusion. To be pre- 
cise, following McKee & Tan [|15J], we shall define Popu- 
lation III stars as those stars with a metallicity sufficiently 
low that it has no effect on either the formation or the 
evolution of the stars. The value of the critical metallic- 
ity for star formation — i.e., the value below which the 



metals do not influence star formation, e.g. the initial 
mass function (IMF) — is uncertain, with estimates rang- 
ing from ^ 1O^^Z0 if the cooling is dominated by small 
dust grains that contain a significant fraction of the met- 
als lUSll to ~ lO^^^Z© if the cooling is dominated by the 
fine structure lines of C and O Ill9ll . It is possible that 
even lower values of the metallicity could significantly 
affect the evolution of primordial stars (Meynet, private 
communication). Among Population III stars, we distin- 
guish between the first and second generations, which 
we term Population III. 1 and III. 2, respectively. We de- 
fine Population III. 1 stars as those for which the initial 
conditions are determined solely by cosmological fluc- 
tuations. By definition, the collapse of gas to form these 
stars is mediated by the cooling properties of primordial 
or near primordial composition gas that has not been sig- 
nificantly influenced by stellar or AGN radiation fields. 
Since we expect the effects of radiation to be more per- 
vasive than the spread and mixing of metals, the compo- 
sition of Pop III.l halos should in fact be precisely that 
resulting from big bang nucleosynthesis. 

As discussed later in this article, the initial conditions 
for Population III. 2 stars are significantly affected by 
other stars, primarily by radiation. Lyman- Werner band 
radiation destroys H2 molecules, reducing the cooling ef- 
ficiency of the gas. Ionizing radiation heats up the gas 
to temperatures ^ 25,000 K ifmUTIl . depending on the 



temperature of the radiation field. While ionized, thermal 
expansion can drive the gas out of low-mass dark mat- 
ter halos. Shocks associated with H II regions and super- 
nova remnants can also ionize and displace the gas. If the 
gas has a chance to recombine, the relatively high resid- 
ual electron fraction catalyzes molecule formation, par- 
ticularly HD formation, which can dramatically enhance 
subsequent cooling, perhaps reducing the characteristic 
star-formation mass 1.20,1 . Greif & Bromm ll2lll use the 
term "Population II.5" to describe stars that form from 
gas in which HD cooling is important, whereas in our 
terminology these would be Population III. 2 stars. We 
believe that it is better to describe all essentially metal- 
free stars as "Population III." It should be noted that this 
definition of Population III. 2 stars includes all Popula- 
tion III stars that were significantly affected by previ- 
ous generations of star formation, even if that did not 
result in significant HD production. According to Greif 
& Bromm |21], Population III.l stars are relatively rare, 
constituting only about 10% by mass of all Population 
III stars. Note, however, that the precise intensity of the 
background radiation field that causes a transition from 
Pop III. 1 to Pop III. 2 remains to be determined from nu- 
merical calculations. 



POPULATION III.l STAR FORMATION 
Initial Conditions and Accretion Rate 

The initial conditions for the formation of the first stars 
are thought to be relatively well understood: they are de- 
termined by the growth of small-scale gravitational in- 
stabilities from cosmological fluctuations in a cold dark 
matter universe. The first stars are expected to form at 
redshifts z ~ 10 — 50 in dark matter "minihalos" of mass 
^ IO^Mq 1I22I1 . In the absence of any elements heavier 
than helium (other than trace amounts of lithium) the 
chemistry and thermodynamics of the gas are very sim- 
ple. Once gas collects in the relatively shallow potential 
wells of the minihalos, cooling is quite weak and is dom- 
inated by the ro-vibrational transitions of trace amounts 
of H2 molecules that cool the gas to ^ 200 K at densi- 
ties «H ~ lO^cm"^ i23ll24ll . As the gas core contracts to 
greater densities, the cooling becomes relatively ineffi- 
cient and the temperature rises to ^ 1000 K. At densities 
^ 10"^ cm^^ rapid 3-body formation of H2 occurs, cre- 
ating a fully molecular region that can cool much more 
efficiently. This region starts to collapse supersonically 
until conditions become optically thick to the line and 
continuum cooling radiation, which occurs at densities 
^ 10'^ cm^^. Recent 3D numerical simulations have ad- 
vanced to densities of order 10^' cm^^ (see contributions 
by Turk et al. and Yoshida et al., these proceedings). 



but have trouble advancing further given the short sim- 
ulations timesteps required to resolve the dynamics of 
the high density gas of the protostar Further numerical 
progress can be achieved by introducing sink particles 
fell. 

We have adopted the alternative approach of modeling 
the protostar's accretion and evolution analytically |26]. 
The accretion rate depends on the density structure and 
infall velocity of the gas core at the point when the star 
starts to form. Omukai & Nishi 12.70 and Ripamonti et al. 
[28] showed that the accreting gas is isentropic with an 
adiabatic index 7~ 1.1 due to H2 cooling; i.e., each mass 
element satisfies the relation P = Kp^ with the "entropy 
parameter" K = const. In hydrostatic equilibrium — and 
therefore in a subsonic contraction — such a gas has a 
density profile p °^ r^^P with kp ~ 2.2. This is confirmed 
by numerical simulations, which show that just before 
protostar formation the gas at the center of the minihalo 
has organized itself into an approximately singular poly- 
tropic sphere (with a modest degree of flattening about a 
rotation axis) with kp — 2.2. In our model, we describe 
the normalization of the density structure of core via the 
entropy parameter 



K' = 



P/P' 



T' 
1.88xlO'2cgs 300 K 



lO'^cm- 



0.1 



(1) 



where Tj^ = T + jj.a^^^y/k is an effective temperature that 
includes the modest effect of subsonic turbulent motions 
that are seen in numerical simulations [i23il . 

Simulations show the gas is inflowing subsonically 
at about a third of the sound speed [23]. We thus use 
Hunter's 129] solution for mildly subsonic inflow (Mach 
number =0.295), with a density that is 1.189 times 
greater than a singular isothermal sphere at f = 0, and 
the accretion rate is 2.6 times greater than the Shu ll30ll 
solution. 

Feedback from the star, whether due to winds, pho- 
toionization, or radiation pressure, can reduce the accre- 
tion rate of the star We define a hypothetical star-|-disk 
mass, nifii Q, and accretion rate, m*^/,o, in the absence 
of feedback. In this case, the star+disk mass equals the 
mass of the core out of which it was formed, m^^/.o = 
M{r). The instantaneous and mean star formation effi- 
ciencies are 



^^d 



C:^d 



Mifd _ m*d 



'w*j,o 



(2) 



(3) 



Assuming the Hunter solution applies for a singular 
polytropic sphere with 7 = 1 . 1, the accretion rate is then 



m^d 



-- O.O26e«//:"^/^(M/M0)-3/^ M© yr-^ (4) 



with the stellar mass smaller than the initial enclosed 
core mass via m* = m^,d/{l +fd)= £*dM/{l +fd)- We 
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FIGURE 1. Mass accretion rate onto the protostar+disk as 
a function of their total mass m»^ for the case of negligible 
stellar feedback. Solid line is the fiducial model from Tan & 
McKee [26] (with K' = 1) from eq. Q. Dotted line is from the 
ID model of Omukai & Nishi |27]. Dashed line is the analytic 
result from Ripamonti et al. 1 28]. Dot-dashed line is the settling 
inflow rate at the final stage of the simulation of Abel et al. |23l, 
now as a function of the enclosed mass. Note that the decline of 
this rate at small masses is due to the lack of the full set of high 
density cooling processes in their simulation. Long-dashed line 
is the equivalent quantity from Yoshida et al. 13111 . Dot-long- 
dashed line is the sink particle accretion rate of Bromm & Loeb 
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choose a fiducial value of /j — 1/3 appropriate for disk 
masses limited by enhanced viscosity due to self-gravity. 
We compare this analytic accretion rate (for the case of 
no feedback) with numerical estimates in Fig. 1. 

O'Shea & Norman |,32,] studied the properties of Pop 
III.l pre-stellar cores as a function of redshift. They 
found that cores at higher redshift are hotter in their 
outer regions, have higher free electron fractions and 
so form larger amounts of H2 (via H^), although these 
are always small fractions of the total mass. As the 
centers of the cores contract above the critical density of 
10'* cm^^, those with higher H2 fractions are able to cool 
more effectively and thus maintain lower temperatures 
to the point of protostar formation. The protostar thus 
accretes from lower-temperature gas and the accretion 
rates, proportional to c^ °^ T^'^, are smaller. Measuring 
infall rates at the time of protostar formation at the scale 
of M = IOOM0, O'Shea & Norman find accretion rates of 
- IO^'^Mq yr-^ at z = 30, rising to - 2 x IQ-^Mq yr"' 
at z = 20. This corresponds to a range in K' of 0.20 
to 2.4. However, for Hunter's mildly subsonic solution, 
which we have suggested is closest to the numerical 
simulations, the accretion rate increases from 0.70c^/G 
at large radii (r 3> c,f, where f = is the moment of 



protostar formation) to 2.58c^/G at small radii (Hunter 
1977), an increase of a factor 3.7. This demonstrates that 
caution should be exercised in inferring accretion rates at 
late times from those measured at early times. 

Spolyar et al. |33] considered the effect of heating 
due to WIMP dark matter annhilation on the collapse 
of Pop III.l pre-stellar cores. In their fiducial model of 
the dark matter density profile and for a WIMP mass 
of 100 GeV, heating dominates cooling (as evaluated 
by Yoshida et al. uM) for nn > lO'"* cm"^, i.e. in the 
inner ~ 20 AU. In this case the collapse can be expected 
to be halted inside this region and a quasi-hydrostatic 
object created, although detailed models of the dynamics 
remain to be worked out. However, as more baryons 
from the outer regions cool and join this central core, 
the luminosity needed to support it increases rapidly 
with mass; the increase provided by WIMP annihilation 
is uncertain. However, in view of the fact that most 
of the later accretion of baryons should occur through 
a disk (see below), which would tend to enhance the 
baryon to dark matter mass ratio in the central disk and 
star, we expect any heating due to WIMP annhilation 
will diminish in importance as the stellar mass grows. 
It should be noted that the WIMP annhilation heating 
rate depends sensitively on the density profile of the 
dark matter, which is so far not adequately resolved in 
numerical simulations of Pop III. 1 core formation. 



Protostar and Accretion Disk Evolution 

Assuming dark matter annhilation does not have a sig- 
nificant effect on the main phase of protostellar accre- 
tion and noting that there may be a range of accretion 
rates depending on the formation redshift |32], we now 
consider what happens to the collapsing gas as it forms 
a protostar and accretion disk at rates that are initially 
;> IQ-^Mq yr-' (i.e. eq. |4]with /:' - 1). 

Rotation of the infalling gas has a dramatic effect on 
the evolution of the protostar, since it leads to lower gas 
densities and optical depths in the regions near the stellar 
surface that are above and below the disk. This leads 
smaller photospheric radii and thus higher temperature 
radiation fields that can then have a stronger dynamical 
influence on the infall. Following the treatment of Tan & 
McKee Ii261 . we parameterize the rotation in terms of 



/Kep 



Vrot(''sp) 
V'Kep('"sp) 



(5) 



the ratio of the rotational velocity to the Keplerian veloc- 
ity measured at the sonic point at r^p. Averaging in spher- 
ical shells, Abel, Bryan, & Norman [23] found /xep '^ 0.5 
independent of radius, so we adopt this as a fiducial 
value. If angular momentum is conserved inside the sonic 



point, then the accreting gas forms a disk with an outer 
radius 



r-d 



1280 ("^y {'!^Y\'-^°n AV, 



1850 



2 9/7 
_<2_ 

0.5 / a:'io/7 



/Kep \ 



AU, 



(6) 



where the -^ is for the case with f^ — 1/3. 

The first gas to collapse has a small disk- 
circularization radius and falls directly on to the 
protostar. The size of the protostar depends on the 
accretion rate during its formation history. At lower 
masses there is a balance in the size set by the need to 
radiate the luminosity, which is mostly due to accretion, 
with a photospheric temperature that is likely to be close 
to ^ 6000 K because the opacity due to H^ rapidly 
declines below this temperature. Under the assumption 
of spherical accretion. Stabler, Palla, & Salpeter Il34ll 
found the protostellar radius to be 



r* ~ 67 (m*/M(: 



nO.27 -0.41 

) '«*.-2 



R, 



0; 



(7) 



where m^.-i = m*/{iO^ Mq yr^ ). For the accretion 
rates typical of primordial star formation we see that the 
size is large, but small compared to the disk radius for 
masses ^Mj,. 

As the star grows in mass its internal luminosity in- 
creases rapidly and begins to dominate over accretion lu- 
minosity [35, 26]. Once the star is older than its instanta- 
neous Kelvin-Helmholz time (i.e. its gravitational energy 
divided by its internal luminosity), it begins to contract 
towards the main sequence configuration. At this time 
there is an associated temporary swelling of the outer 
layers of star by a factor of about two due to structural re- 
arrangements inside the star. According to Schaerer Ii36i1 . 
the main sequence radius for nonrotating stars is 



r*~4.3m^;2 ^© (main sequence) 



(8) 



to within 6% for 0.4 < m*.2 < 3. If the star continues to 
accrete and gain mass, it will approximately follow this 
mass-size relation, although rotation will tend to swell 
the equatorial surface layers somewhat. The evolution 
of the protostellar radius is shown for several models in 
Fig.E] 

The high accretion rates of primordial protostars make 
it likely that the disk will build itself up to a mass that 
is significant compared to the stellar mass. At this point 
the disk becomes susceptible to global (m = 1 mode) 
gravitational instabilities i38ll39ll . which are expected to 
be efficient at driving inflow to the star, thus regulating 
the disk mass. Thus we assume a fixed ratio of disk to 
stellar mass, fd — 1/3, in our models. 

Accretion through the disk may also be driven by lo- 
cal instabilities, the effects of which can be approximated 
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FIGURE 2. Evolution of protostellar radius with mass, based 
on the model of Tan & McKee 12611 . All models shown have 
a rotation parameter /Kep = 0.5. The solid line shows the 
evolution of the fiducial model with A"' = 1, a Shakura-Sunyaev 
disk viscosity parameter of Ofss = 0.01, and a reduction of 
accretion rate due to feedback effects |15], which becomes 
important for ;n, > 50Mr;). The star joins the main sequence 
1 36], shown by the /ong dashed line, at about SOM©. The dot- 
dashed line, visible only from 50— lOOM,;,, shows the same 
model but with no reduction in accretion rate due to feedback 
effects. The dot-long-dashed line, visible up to about 2OM0 
shows the fiducial model, but with 0!ss = 0.3 (appropriate for 
viscosity driven by self-gravity |37]). The dotted line shows the 
fiducial model but with K' — 2, while the dashed line shows the 
K' = 0.5 case. 



by simple Shakura-Sunyaev ttss-disk models. Two di- 
mensional simulations of clumpy, self-gravitating disks 
show self -regulation with ajs ~ (ilfth)"' up to a maxi- 
mum value ttss — 0.3 113711 . where Q. is the orbital angular 
velocity, ftj, = T.kTi.^/{aT^ff^) is the thermal timescale, 
T, is the surface density, T^d is the disk's central (mid- 
plane) temperature, and Tiff.d the effective photospheric 
temperature at the disk's surface. Fragmentation occurs 
when ilfth < 3: this condition has the best chance of be- 
ing satisfied in the outermost parts of the disk that are still 
optically thick. However, Tan & Blackman 1,40] consid- 
ered the gravitational stability of constant ttss = 0.3 disks 
fed at accretion rates given by eq. |4] and found that the 
optically thick parts of the disk remained Toomre stable 
(Q> I) during all stages of the growth of the protostar. 
We therefore expect that most Pop III. 1 accretion disks, 
at least in their early stages, will grow in mass and mass 
surface density to the point at which gravitational insta- 
bilities, both global and local, then mediate accretion to 
the star 

In addition to gravitational instabilities, the magneto- 
rotational instability (MRI) may develop if dynamically 



important magnetic fields are generated by dynamo ac- 
tion in the disk (see below). This process could be partic- 
ularly important in the inner accretion disk, where grav- 
itational instability is suppressed. The MRI is thought 
to yield viscous stresses that correspond to somewhat 
smaller values of ass ^ 0.01 |i4li1. although the actual ef- 
fective value remains unclear. 



Generation of Magnetic Fields in the Accretion Disk 

Given the right conditions, dynamo action can am- 
plify a seed magnetic field. The basic requirement is that 
the field is frozen into the gas and that gas motions re- 
peatedly stretch the field, e.g. via turbulent or shearing 
motions. Here we briefly describe the results of Tan & 
Blackman |!40!], who considered these issues. 

The seed field for the primordial star-forming environ- 
ment can be estimated from the Biermann battery mech- 
anism for non-barotropic flows. On scales much larger 
than that of an individual star, numerical simulations in- 
dicate that by z = 18, seed fields of order Bq ^ 10^^^ 
Gauss can be generated |42]. We assume such a field is 
generated inside the region around the minihalo on scales 
^ lOOpc. The infall speeds within this region are of order 
1 km s^', vastly greater than the drift speeds of the ions 
through the neutrals (the ionization fractions are much 
greater than the equilibrium values for the gas, which has 
a minimum temperature of ^ 200 K). Thus the field is 
frozen into the contracting gas: going from parsec to AU 
scales increases the field strength by a factor of '^ lO'", 
so the seed field in the disk may then be ~ 10^'^G. This 
estimate is highly uncertain. Seed fields several orders of 
magnitude larger could arise in the protostellar disk from 
the Biermann battery term if the pressure and density 
gradients are misaligned inside this region as well. A fur- 
ther complication that could lead to lower field strengths 
is the possibility of efficient turbulent diffusion of the 
field so that flux freezing and field advection do not play 
as significant a role. However, since the disk-dynamo 
amplifies the initial field exponentially, the uncertainties 
in the strength of the seed field do not significantly affect 
the time to reach saturation. 

A seed field is expected to be amplified by a combi- 
nation of shearing and turbulent motions in the disk. As 
discussed above, initially the turbulence is thought to be 
due to self-gravity |37], but as the field strength grows, 
this could be replaced by the MRI. The field strength 
grows exponentially with a growth time ^ ttss Qr^ . 
Field growth stops when the fields become dynamically 
important, i.e. forB ^ (47i:p)''^assCi. The generation of 
strong B-fields that are ordered on scales large compared 
to the disk thickness requires turbulence in a stratified 
disk: the helical disk dynamo[43]. More specifically the 



turbulence must generate net helicity in each hemisphere. 
If this occurs, and numerical simulations are required for 
confirmation, then these large-scale B-fields could start 
to drive a hydromagnetic bipolar disk wind, whose feed- 
back effects will be discussed below. 

The presence of strong B-fields in the disk and their 
possible advection into the protostar, then allows a star- 
disk interaction that can help the star to shed angular 
momentum and may set the rotation rate with which 
it is born. Note that models of primordial protostellar 
evolution predict that the star is mostly radiatively stable 
in the pre-H nuclear burning phase, which implies there 
is very little dynamo amplification of B-field in the star 
independent of the disk. 



Radial Structure of the Accretion Disk 

Given these considerations, a first approximation for 
the radial structure of the protostellar accretion disk can 
be made by using the standard theory of steady, thin, vis- 
cous accretion disks, with a spatially constant viscosity 
parameter, ttss, and for simplicity ignoring energy injec- 
tion from the star. Tan & McKee |26] calculated the evo- 
lution of the radial structure, assuming the disk is fed at 
a rate given by equation (HJi and the inner boundary is at 
the protostellar surface, r*. The viscosity is assumed to 
be a function of the total pressure. Typically gas pressure 
dominates over radiation pressure in the earlier stages 
(m* < 2OM0). We chose ttss = 0.01 as a fiducial value, 
that may arise from viscosity provided by the magneto 
rotational instability IHTT . 

The radial structure of the disk is governed by the 
equations of energy conservation and of angular momen- 
tum conservation. Energy conservation gives the emer- 
gent flux as 
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Here 
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/=1 



1/2 



(9) 



(10) 



(11) 



is the factor that embodies the boundary condition that 
angular momentum cannot be transferred across a sur- 
face on which the angular velocity has no gradient; CJo is 
the cylindrical radius at which dH/dU} vanishes, which 
we take to be equal to the stellar radius. The dimension- 
less factor 0/ describes the advection of thermal and in- 
ternal energy in the disk and is generally less than unity. 
To evaluate the angular momentum transfer in the disk, 
we adopt the a-disk model of Shakura & Sunyaev l44ll . 



in which the transverse stress in the disk is proportional 
to the pressure, wojii, = —j(^P (we have included the 
factor I to conform with the convention of Frank et 
al. 1451 ). The equation describing angular momentum 
transport is then 



(12) 



m*i2/ = 6Ka / Pdz, 



where Zs is the height of the surface of the disk. 

The radial disk structure at three different stages (m* = 
1, 10, lOOMg) is shown in Fig. [51 with the radial coordi- 
nate measured in units of the protostellar radius. 



Vertical Structure of the Accretion Disk 



A numerical solution of these equations using the 
opacities of Iglesias & Rogers ||48I1 is shown in Fig. |4] 
for one particular stage of the protostellar evolution for 
a location at 10/?* around a lOOM© main sequence star, 
accreting at 2.4 x lO^^MQyr^', i.e. the fiducial rate from 
a /T' = 1 core with no reduction due to feedback. 

More generally, we follow the disk structure during the 
course of the protostellar evolution (i.e. as m*, m* and r* 
evolve). As we discuss below, the disk thickness is one 
important factor in determining the stellar mass at which 
accretion is shut off. At masses m* ^ lOOMp we find 
the aspect ratio of the disk surface, Zs/r, is about twice 
h/r, where h is the scale height, and has a very weak 
dependence on radius. 



McKee & Tan 11511 calculated the vertical structure 
of Pop III.l accretion disks. The vertical structure of 
the disk is governed by three equations: first is the first 
moment of the radiative transfer equation, 

dPr^i 



pKpF 



dz 



(13) 



where Kp is the flux-weighted mean opacity per unit 
mass and F{z) is the radiative flux. We assume that 
the effective optical depth for true absorption, T* = 
(fabsfscatt)''^, IS much greater than unity so that the gas 
is approximately in LTE [44, 46]. Then Pjad — ^aT'^ and 
Kf ~ (Cr, where jcr is the Rosseland mean opacity per 
unit mass, so that equation ( fT3] l reduces to the equation 
of radiative diffusion. 






3krpF 



(14) 



The second equation describes the growth of the flux 
due to viscous dissipation |47], 



dF 



do. 9 

4^ 



— = -(piwnj^dj — = -<^iaQ.P. 



(15) 



We have included the factor 0/ to allow for the reduction 
in the flux by the advection of internal energy. In addition 
to the factor j 0/, equation (flST l differs from the expres- 
sion adopted by Shakura & Sunyaev 04411 in that it has 
dF /dz °^ P rather than <=<: p . One can show, however, that 
the height of the disk is very insensitive to this change. 
Integration of equation ( fTSl ) together with equation ( fT2l i 
leads directly to the energy equation ( fTOl i. 

Finally, we have the equation of hydrostatic equilib- 
rium, 

dp _ pGnifZ 

where the pressure P includes both gas pressure and 
radiation pressure. 



(16) 



P = P,+Pmd = 



pkT Aar 



M 



3c 



(17) 



Feedback Processes and the IMF 

Radiative Feedback 

The evolution of the bolometric luminosity of the pro- 
tostar and the inner accretion disk is shown in Fig. |5] 
for three different accretion rates corresponding to K' = 
0.5, 1,2. The luminosities are sub-Eddington for almost 
the entire evolution. Only the K' = 2 model (with the 
highest accretion rate) exceeds the Eddington limit: at the 
low-mass limit this is due to the initial condition having a 
stellar radius that is too small and the star quickly adjusts 
by undergoing rapid expansion; at about SOM© the star is 
close to joining the main sequence and so has a relatively 
small size and is still accreting relatively efficiently be- 
cause its H II region is still confined (see below), thus 
resulting in a high accretion luminosity. However, in this 
latter case we do not expect accretion to be significantly 
affected by radiation pressure because much of the lumi- 
nosity will be directed in polar directions and most of the 
accretion is of neutral gas from the equatorial regions. 

As discussed in McKee & Tan lUSll . we expect ioniz- 
ing feedback to be more important than continuum radia- 
tion pressure at reducing accretion efficiency and eventu- 
ally setting the final stellar mass. The H-ionizing photon 
luminosities for the K' — 0.5, 1,2 models are shown in 
Fig. |6l At lower masses, the protostars are much larger 
than the corresponding main sequence stars and so have 
cooler photospheric temperatures and smaller ionizing 
luminosities. Once the protostar is older than about its 
instantaneous Kelvin-Helmholz time, which occurs at 
higher masses for higher accretion rates, the star con- 
tracts towards the main sequence and the ionizing lumi- 
nosity rises accordingly. 

Figure |7] gives an overview of the ionizing feedback 
processes occurring near the protostar At early stages 
the ionizing flux from the protostar is smaller than the 
flux of accreting neutral atoms to its surface and the H II 
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FIGURE 3. Protostellar disk structure 1 40] for models with ttss = 0.01 and m* = 1 , 10, IOOMq , for which r, = 100, 300, 4Rq and 
m* = (17,6.4,2.4) x IO^^Mq yr^', respectively, i.e. from eq.|4]with no feedback. From top to bottom the panels show (1) surface 
density, Z; (2) ratio of scaleheight to radius, h/r, and ratio of gas pressure to total pressure, j3; (3) midplane ionization fractions 
of H+, He+, He +; (4) disk midplane temperature, T^j^ (the dotted lines show results for when the ionization energy is neglected); 
(5) number of orbits before accretion to the star, n^j.^,, Toomre Q stability parameter, and Rosseland mean opacity k, evaluated at 
the midplane. Note that all quantities are azimuthal and temporal averages of the disk, which, being turbulent, would exhibit local 
fluctuations. 
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FIGURE 4. Vertical structure of the accretion disk at r = 
lOr* ~ 43Rq for m, = 100 Mq, K' = 1, /xep = 0.5 and no 
reduction in accretion efficiency due to feedback. 
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FIGURE 5. Protostellar bolometric luminosities for models 
with A"' = 0.5, 1,2 (dashed, solid, dotted lines). In each case the 
total luminosity is shown with the higher line and the accretion 
luminosity (boundary layer + inner disk) is shown with the 
lower line. The Eddington limit is indicated with the dot-lone- 
dashed line and the zero age main sequence luminosity I3al 
with the long-dashed line. 



region is confined there. As the ionizing luminosity in- 
creases, the H II region can begin to expand into the in- 
falling envelope, penetrating furthest along the rotation 
axes. At the ionized-neutral boundary, radiation pressure 
feedback is exerted due to resonant scattering of FUV 
radiation in the Lyman-a damping wings. As a result 
of the high column densities of neutral gas around the 
H II region, this radiation is trapped and the pressure am- 



FIGURE 6. Protostellar H-ionizing photon luminosities for 
models with K' = 0.5,1,2 (dashed, solid, dotted lines). The 
total luminosity is shown with the higher line and the accretion 
luminosity (boundary layer + inner disk) is shown with the 
lower line. The zero age main sequence H-ionizing luminosities 
I3al are indicated with the long-dashed line. 



plified by large factors. For typical rotation parameters 
/Kep — 0.5, this radiation pressure becomes larger than 
the ram pressure of the infalling gas in the polar direc- 
tions for stellar masses of order 2QMq. However, once 
the infall is reversed at the poles, the Lyman-a photons 
can escape and the accretion in other directions proceeds 
relatively unimpeded. Thus we expect Lyman-a radia- 
tion pressure to have only a minor effect on the accretion 
efficiency. 

The next feedback effect to occur is the expansion of 
the H II region to distances larger than the gravitational 
escape radius r^ for ionized gas. This distance is 



GipEddm- 



*d 



2.5 



260fed — '««/,2 AU, (18) 



where T4 is the ionized gas temperature in units of lO'* K 
and we have taken the gravitating mass to be m*^ and 
we have allowed for the decrease in the radius due to 
radiation pressure from electron scattering through the 
factor 

(/)Edd = 1 - 7^ , (19) 

i-Edd 

where Lgdd = 47rGOTc/K'xhomson is the Eddington limit. 

Shortly after the H II region reaches r^, pressure forces 
in the ionized gas become large enough to reverse the 
infall to the protostar This effect will occur first in the 
polar directions, typically at about 5QMq in our fiducial 
model. For rotation parameters /Kep — 0.5, expansion in 
the equatorial directions (just above the accretion disk) 






HII region 
breakout 




FIGURE 7. Overview of the accretion geometry and feedback processes involved in primordial star formation, (a) Top left: Cross 
section of the accretion geometry: the dashed lines show streamlines of the rotating, infalling gas, with figure of revolution from 
each streamline separating 10% of the total infall from this hemisphere. The aspect ratio of the accretion disk is realistic, while the 
size of the star has been exaggerated for clarity, (b) Top right: The red shaded region shows the extent of the H II region, which 
at this relatively early stage is still confined inside the gravitational radius for the escape of ionized gas, rg. Lyman-a radiation 
pressure feedback should be strong enough to prevent accretion in the polar directions, (c) Bottom left: The stellar mass and 
ionizing luminosity have grown, and the H II region is just in the process of breaking out of the accretion flow. Once a significant 
volume beyond rg is ionized, accretion from these directions is expected to be shut off. (d) Bottom right: Final stage of accretion 
involves shadowing of the equatorial region by the disk, which at the same time is photoevaporated. The competition between this 
photoevaporative outflow and the residual accretion rate sets the final mass of the star. 
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FIGURE 8. Mass scales of H II region break out as a 
function of the rotation parameter /Kep- The lower dashed line 
marked "Polar" shows the mass scale of the protostar at which 
the H II region reaches rg (based on star plus disk mass) along 
the rotation axis of the protostar The upper dashed line marked 
"Equatorial" shows the mass scale of the protostar when the 
H II region reaches rg in a direction just above the disk plane 
(0.9;r/2 from the rotation axis). 



occurs very soon after polar breakout (see Fig. [8] and 
Table [T]i. Higher accretion rates and smaller rotation 
parameters result in larger mass scales for H II region 
breakout. 

The protostar and accretion envelope should evolve to 
a state in which the entire region above and below the ac- 
cretion disk has been ionized (Fig.|7}l). No accretion will 
occur from these directions. The accretion disk, which is 
still mostly neutral in its outer regions, has a finite thick- 
ness and will shield an extended equatorial region from 
direct ionizing flux from the protostar, allowing accretion 
to be able to continue from these regions. The reduction 
in accretion rate caused by the expansion of the H II re- 
gion and the fraction allowed in the disk-shadowed re- 
gion is shown in Fig. |9] where the fiducial no-feedback 
case is also shown. 

Also shown in Fig. |9]is the mass loss rate due to the 
process of disk photoevaporation ||49I1 . Ionization from 
the protostar creates an ionized atmosphere above the 
neutral accretion disk, which then scatters some ionizing 
photons down on to the shielded region of the outer 
disk, beyond r^. An ionized outflow is driven from these 
regions at a rate 



*evap 



:4.1 X 10 



-5r.l/27,0.4^1/2 



J 49 -14 



d.2 



Mgyr- 



(20) 



where 549 is the H-ionizing photon luminosity in units of 
10^^ photons s^^ When this rate becomes comparable 
to the accretion rate to the disk, then we expect further 



FIGURE 9. Feedback-limited accretion: fiducial case. The 
evolution of the accretion rate versus protostellar mass is shown 
for the fiducial model (/Kep = 0.5, K' = l,Ti = 25,000 K) in 
the cases of "no feedback" and "with feedback". In the latter, 
the accretion efficiency is reduced as the H II region expands to 
rg and beyond. However, accretion is allowed to continue from 
directions that are shadowed by the disk photosphere. The disk 
structure and protostellar structure and feedback are calculated 
self-consistently given the evolution in lii:^. Also shown is the 
photoevaporative mass loss rate, mevap, which starts once the 
H II region has broken out in the equatorial direction and grows 
as the ionizing flux increases. We see that this mass loss rate 
becomes greater than the accretion rate at m^, ~ 137M©, and 
we identify this mass scale as our best estimate of initial mass 
scale of the first stars. 



growth of the stellar mass to be very limited. For sim- 
plicity we equate the final stellar mass to that mass at 
which the photoevaporative mass loss rate equals the ac- 
cretion rate. It is about 14OM0 in the fiducial case, and 
Table [T] and Fig. [TO] summarize other cases. 

An analytic estimate can be made for the final proto- 
stellar mass. Assuming the H-ionizing photon luminosity 
is mostly due to the main sequence luminosity of the star 
we have 



: 7.9 X 10'*'^ ^s'nli 



ph s 



-1 



(21) 



which for ^5 = 1 is a fit to Schaerer's (2002) results 
for the ionizing luminosity of main sequence primordial 
stars; the fit is accurate to within about 5% for 6QMq <^ 
nit, <^ 3OOM0. Then the photoevaporation rate becomes 



^evap 



1.70xlO-Vy'(l+/rf)'/' 



T, X 0.4 

h \ 5/4 
15 



Meyr 



-1 



(22) 



The accretion rate onto the star-disk system is given by 
equation (|4|i. Equating this with equation ( |22] |. we find 



TABLE 1. Mass Scales of Population III. 1 Protostellar Feedback 

K' /Kep TiliW^K) m.^iMat. ™*,eb(M0ii m*,evap (MoJ^ 



1 0.5 


2.5 


45.3 


50.4 


137£ 


1 0.75 


2.5 


37 


41 


137 


1 0.25 


2.5 


68 


81 


143 


1 0.125 


2.5 


106 


170 


173 


1 0.0626 


2.5 


182 


330£ 


256 


1 0.5 


5.0 


35 


38 


120 


1 0.25 


5.0 


53.0 


61 


125 



0.5 


0.5 


2.5 


23.0 


24.5 


57 


2.0 


0.5 


2.5 


85 


87 


321 



* Mass scale of HII region polar breakout. 

t Mass scale of HII region near-equatorial breakout. 

** Mass scale of disk photoevaporation limited accretion. 

* Fiducial model. 

^ This mass is greater than hi, evap in this case because it is calculated without allowing 
for a reduction in m, during the evolution due to polar HII region breakout. 



that the resulting maximum stellar mass is 



Maxm. 



*/:2 



6.3 



eT-eT<^TK'^"'' 



{l+fd)^'"' 



2.5 



0.24 



(23) 

Recall that g^.^/ is the instantaneous star formation 
efficiency — i.e., the ratio of the accretion rate onto the 
star to the rate that would have occurred in the absence 
of feedback. Here this ratio is just the shadowing factor, 
/sh, i.e. the fraction of the sky as seen from the protostar 
that is blocked by the disk. For this simple analytic es- 
timate, we make the approximation that the shadowing 
sets in when the stellar mass reaches mi, so that £«/ = 1 
until the mass of the central star reaches my and e^,d = /sh 
thereafter It is then straightforward to show that 



C'^d 



./sh 



l-(l-/sh)(wi/w«/)' 



(24) 



provided that m^,d >m\. Note that the average efficiency 
Eifd = 1 at the onset of shadowing (m^^ — m\) and that 
e*c/ ~^ ./sh at late times m^^ ^ 'Wi- Normalizing /sh to a 
typical value of 0.2 we find 



Max m 



*/,2 



= 1.45/^'^''/4^ 

r \ 28/47 
Jsh ^ 

0.2 



2.5 
025 



0.24 



12/47 



(25) 



where we have set the ionizing luminosity factor ^s = 1 
and the disk mass fraction // = i; we have normalized 
e^,d to a value of 0.25, which is approximately correct 
for K' = 1 and for wi ~ 50Mq and m«/ = 2OOM0. This 
analytic estimate therefore also suggests that for the fidu- 
cial case (K' =1) the mass of a Pop III.l star should be 
~ 14OM0. 



The uncertainties in these mass estimates include: (1) 
the assumption that the gas distribution far from the star 
is approximately spherical — in reality it is likely to be 
flattened towards the equatorial plane, thus increasing the 
fraction of gas that is shadowed by the disk and raising 
the final protostellar mass; (2) uncertainties in the disk 
photoevaporation mass loss rate due to corrections to the 
Hollenbach et al. [49] rate from the flow starting inside rg 
and from radiation pressure corrections; (3) uncertainties 
in the H II region breakout mass due to hydrodynamic 
instabilities and 3D geometry effects; (4) uncertainties in 
the accretion rate at late times, where self-similarity may 
break down [25]; (5) the effect of rotation on protostellar 
models, which will lead to cooler equatorial surface tem- 
peratures and thus a reduced ionizing flux in the direction 
of the disk; (6) the simplified condition, ;«evap > »?*(/, 
used to mark the end of accretion; and (7) the possible 
effect of protostellar outflows (discussed below). 



Mechanical Feedback 

Protostellar outflows, thought to be launched by large- 
scale magnetic fields threading the inner accretion disk, 
are ubiquitous from present-day protostars. The momen- 
tum flux in these flows typically exceeds that due to 
radiation pressure by several orders of magnitude. As 
described above. Tan & Blackman 1.40.1 discussed how 
turbulence in a stratified disk that generates helicity is 
the necessary condition for production of dynamically- 
strong large-scale magnetic fields by dynamo action. 
They then considered the effect of such an outflow on 
the surrounding minihalo gas, following the analysis 
of Matzner & McKee IdOII . The force distribution of 
centrifugally-launched hydromagnetic outflows is coUi- 
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FIGURE 10. Feedback-limited accretion: effect of ionized 
gas temperature and accretion rate. The fiducial model (/Kep = 
0.5, K' = 1, Ti4 = 2.5 K) shown in Fig. |9] is compared to 
models in which one parameter has been changed: a model 
with Ti4 = 5 and two models with K' = 0.5,2. The dashed 
lines show the accretion rate to the star, m*, and the solid lines 
show the photoevaporative mass loss rate, raevap- The change 
in temperature causes relatively minor differences, while the 
change in K', equivalent to a change in rii^ of factors of 4.4 
above and below the fiducial level, leads to roughly a factor of 
2.4 change in the final stellar mass. Note the increase in m* 
for the K' = 0.5 case at around 35M(;, is due to a thickening of 
the inner accretion disk as the star contracts down to its main 
sequence configuration and assumes material at large distances 
still remains to be accreted in the enlarged shadowed region. 



speed of the core, Vesc.c- The location of the surface 
of Pop III.l cores is not particularly well defined: the 
fiducial core of Tan & McKee |26], based on the sim- 
ulation results of Abel, Bryan, & Norman |23], has 
about 1000 Mq of virialized baryonic material with a 
few xlO^M© of surrounding gas still infalling. These 
cores have Vesc.c = 3.22(M/lOOOM0)->/^/:'V7 km/s. 
To find the speed of the swept-up gas, the analysis 
assumes thin, radiative shocks, purely radial motion 
and monopole gravity. We consider cores with angular 
mass distributions of the form dM/dO. = (l/4;r)g(^)M 
with Qin) = (1 - ii^)"/Jq{1 - pi^)" dpi, and eval- 
uate models with n, ranging from 0, isotropic, to 4, 
which describes a flattened distribution that mimics the 
effect of some rotational support. In the notation of 
Matzner & McKee OSOll . the escape condition is given 
by (1 + e2 - Ale^,je(Mesc) = m*/(XM), where X = 
O.132cJln(2/0o)/ln2OO]vesc,c.5[(Pw/«i*)/4Okms-i]-i, 
and Cg is a factor of order unity that accounts for the 
effects of gravity on shock propagation. We estimate 
Cg ci 4.6, for steady winds that decouple from the swept- 
up shell at the core edge in a core with kp = 2.2 (which 
is the fiducial value of Tan & McKee Ii26tl ). We ignore 
the influence of material beyond the core "edge", which 
is partly balanced by our neglect of the wind's influence 
on the shell beyond this point. 

The star formation efficiency due to protostellar out- 
flow winds is given by 



C-w 



1 



l + (m„/OT*)(l -^„,) Jo 



Mcsc 



Gd/i, 



(27) 



mated along the rotation axes, but includes a signifi- 
cant wider-angled component. Matzner & McKee \5V\ 
showed that far from the star the angular distribution of 
the momentum in a radial hydromagnetic wind can be 
approximated by 
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4;rln(2/0o)(l 



^0^- 



-M^ 



(26) 



where ji — cos0 and 9o is a small angle, which is esti- 
mated to be ^ 0.01 for winds from low-mass stars. The 
total momentum of the outflow, /j,^, is calculated by as- 
suming a fixed ratio of mass injected into the outflow rel- 
ative to mass accretion to the star, ih^./ihi = 0.084 in the 
fiducial model of Tan & Blackman. The speed of the out- 
flow is set equal to the escape speed from the stellar sur- 
face. As the star grows in mass and contracts to the main 
sequence the cumulative specific momentum of the out- 
flow per stellar mass rises from pn/m^ ~ 7kms^^ when 
m* ~ IOM0 to Piv/nit 2± 100km s"^ when w* ~ IOOMq. 
To evaluate the efficiency of star formation from a 
core we find the angle, 0esc = cos^'/Xesc, where the 
wind sweeps up core material to the surface escape 



where 0„, = fQ^'"'Pdjj. and P(/i) is the dimensionless 
force distribution in angle. Note that this definition of 
efficiency includes the effect of diversion of a fraction 
m,v/m* — 0.084 of the stellar accretion rate into out- 
flow material. Equation [27] can be evaluated as a func- 
tion of m, (see Fig.fTTT). 0esc gradually increases and Ew 
decreases as the star grows and injects more momentum 
into its surroundings. However, note that Ew is an upper 
limit to the star formation efficiency if it is evaluated be- 
fore the final stellar mass is reached since some material 
at > 0esc would later be ejected. Note also that the sim- 
ple models shown here do not self-consistently account 
for the effects that a smaller mass accretion rate would 
have on protostellar evolution and outflow strength. Nev- 
ertheless they show that by the time m^, has reached a 
hundred to a few hundred solar masses, the outflows have 
been able to eject about half the initial core material. 

Figure [TT] also shows the instantaneous accretion ef- 
ficiency to the star+disk, E^ii, due to radiative feedback 
and disk-shadowing for the fiducial K' — I, /xep — 0.5, 
71,4 = 2.5 model. Comparison of these efficiencies is 
complicated by the fact that, as mentioned, for E^ at a 
given m-t, some gas at angles > 0esc will eventually be 



ejected if the outflow continues. Nevertheless, we see 
that radiative feedback and disk-shadowing are more im- 
portant, i.e. lead to smaller efficiencies, than mechanical 
feedback for masses ^ SOMg in the fiducial case. The 
outflow would need to operate until m^, ^ IOOOMq to re- 
duce £n. to the value of £„j/ seen at m* ~ lOOM©. The 
possibility of lateral deflection of gas streamlines by the 
outflow is not treated by the above analysis, but we ex- 
pect this would tend to increased e,,,. This effect can be 
studied with numerical simulations. Machida et al. ll52ll 
have carried out 3D MHD simulations of the very earliest 
stages of protostellar outflow feedback, up to times when 
the protostar has just formed with m* <, O.OlMo. Such 
simulations need to be advanced to later stages to explore 
the effects of a realistic 3D geometry on the outflow-core 
interaction and star formation efficiency. 

While mechanical feedback is unlikely to dominate 
over radiative feedback, it may influence details of the 
radiative feedback processes. Even a small outflow cav- 
ity when OT* ~ 20M(;) will help dissipate Ly-a photons 
and prevent a build up of dynamically important radia- 
tion pressure at this stage. An outflow will tend to facili- 
tate H II region breakout in polar directions by lowering 
gas densities and injecting additional momentum. As it 
is launched from the surface of the disk, a dense bipo- 
lar outflow can confine the protostellar ionizing flux in 
equatorial directions and prevent it reaching outer disk to 
cause photoevaporation [531. However, Tan & Blackman 
140] showed that the fiducial outflow would not be dense 
enough to significantly shield the disk from the ionizing 
luminosities of primordial protostars once m* ^ AQMq. 
On the other hand, the outflow is likely to impact the 
structure of the ionized disk atmosphere and thus affect 
the photoevaporation mass loss rate, perhaps similar to 
the "strong stellar wind" case of HoUenbach et al. 1.49.] . 



Impact of Pop III.l Stars for Larger- Scale 
Feedback and Metal Production 

The predicted final protostellar masses (i.e. the initial 
stellar masses) for the range of model parameters we 
have considered are summarized in Table [T] The stars 
are all "massive" (i.e. m* > SM©) and will thus have sig- 
nificant radiative and mechanical feedback on their sur- 
roundings, and the potential for substantial metal enrich- 
ment via core-collapse or pair-instability supernovae. 

Ionization, outflow, wind and supernova feedback can 
lead to the disruption and unbinding of gas from nearby 
minihalos In particular ionization is likely to be the first 
and most important feedback effect to initiate unbinding: 
ionizing radiation heats up the gas to temperatures ^ 
25,000 K Im UTIl . depending on the temperature of 
the radiation field, with corresponding sound speeds that 
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FIGURE 11. Evolution of star formation efficiency, £„,, due 
to erosion of a lOOOMp gas core by protostellar outflows winds 
for the fiducial /Kep = 0.5 and K' =\ case, with Ofss = 0.01, and 
ignoring radiative feedback processes |40]. The density profile 
of the initial core is specified by dM/dO. = {i/47t)Q{il)M, 
with /I = COS0 and Q{^) = (\ - ^^)" / Jo{l - ^^)" d^. Solid 
line is n = (isotropic core), dotted is « = I, dashed is « = 2, 
long-dashed is n = 3, dot-dashed is « = 4. The dot-long-dashed 
line shows the instantaneous efficiency, e^j, due to radiative 
feedback processes (i.e. H II region breakout, limited by disk- 
shadowing) for the fiducial model (/Kep = 0.5, K' = I, Tn = 
2.5), indicating their greater importance at m* > 50Mq. 



are much larger than the few km/s typical of minihalos. 
Note, however, that in addition to unbinding gas from 
minihalos, the pressure forces associated with ionized 
regions can shock-compress adjacent neutral regions and 
promote gravitational collapse, and this possibility is 
discussed in the next section. Stars with masses in the 
range shown in Table [T] are able to completely ionize 
at least their local minihalo, and more typically a much 
large volume of the surrounding intergalactic medium. 
For recent numerical simulations of photoevaporation of 
minihalos see Whalen et al. [54] and references therein. 
Lyman-Werner band radiation destroys H2 molecules, 
reducing the cooling efficiency of the gas. We can es- 
timate the distance over which star formation is sup- 
pressed from the work of Glover & Brand 155]. Assum- 
ing that the core is in approximate hydrostatic equilib- 
rium and is characterized by an entropy parameter K, we 
find [15] that the time to dissociate H2 is less than the 
free-fall time % if the core is within a distance 
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f Siy/ 10 ^ /abs/diss 
1^1049 s-1 X2 0.01 
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^21/40^,1/4 



PC, 



(28) 

of the protostar, where S'lw is the photon luminosity in 
the Lyman-Werner bands, X2 is the fractional abundance 



of H2, /abs is the fraction of the Lyman-Werner flux 
absorbed by the H2, /diss is the fraction of absorptions 
that resuh in dissociation, and h is the mean density of 
H nuclei. Thus, even a lOOMo star, which has S'lw — 
10"*^^ s^', can suppress star formation in an existing core 
only if the core is relatively nearby. A more detailed 
analysis by Susa Ipq] co mes to similar conclusions (see 
also Ahn & Shapiro [57] and Whalen et al. i54tl ). 

If gas ionized by the Pop III. 1 star has a chance to re- 
combine, e.g. once the star completes its nuclear burning, 
the relatively high residual electron fraction catalyzes 
molecule formation, particularly of HD, which can dra- 
matically enhance subsequent cooling ||20|] . This is the 
most likely scenario for the next generation of star for- 
mation. Population III. 2 (discussed in the next section). 

Lacking metal lines in their atmospheres, main- 
sequence and post-main-sequence Pop III stars are 
expected to have quite weak stellar winds IdSII . This also 
implies their mass at the end of stellar evolution should 
be similar to their initial mass. Rotation is expected to 
enhance mass loss 115911 . although the final mass is still 
thought to be a significant fraction of the initial mass. 

Thus the range of final protostellar masses predicted 
by models of radiative and mechanical feedback (Ta- 
ble [T]i imply a range of final pre-core-coUapse stel- 
lar masses that potentially overlap the range of masses 
necessary to produce pair instability supernovae, 140 — 
260 M0 in the models of Heger & Woosley f6^. Rota- 
tion may lower these limits (S. Woosley, private comm.). 
The lack of the expected nucleosynthetic signature of 
such supernovae in the abundance patterns of very metal 
poor halo stars 1 6 1 ] , could indicate that such massive Pop 
III. 1 stars were relatively rare or that they tended to en- 
rich regions not probed by typical halo stars, perhaps the 
centers of larger galactic halos. The conclusion by Scan- 
napieco et al. |62] that Pop III star formation should be 
fairly widespread in regions now probed by Galactic halo 
stars, can be reconciled with the abundance pattern ob- 
servations if most of this star formation leads to either 
Pop III.l stars from relatively low entropy (K' <^\) gas 
cores or Pop III. 2 stars (discussed below) that also have 
a mass scale below the pair instability threshold (see also 
the study by Greif & Bromm |21]). Further work is re- 
quired to determine the range of pre-stellar core param- 
eters, primarily K' and /xep, exhibited in cosmological 
simulations, in order to predict the frequency of pair in- 
stability supernovae. 



POPULATION III.2 STAR FORMATION 

Population III. 2 stars are those with near primordial com- 
position, so that their formation and stellar evolution are 
independent of metallicity, but with initial conditions 
"significantly" affected by radiative and/or mechanical 



feedback from previous stellar generations (that could be 
Pop III.l, Pop III.2, or even Pop II or Pop I) or AGN. 

We expect radiative feedback (especially FUV de- 
struction of H2) to be far more pervasive than the direct 
mechanical feedback associated with protostellar out- 
flows, winds or supernovae, although shocks driven by 
D-type ionization fronts will be even more penetrating 
than ionization into dense regions. Such shock compres- 
sion may be important in promoting star formation in 
neighboring minihalos of intermediate central density, 
lOcm^^ ^ "H < lOOOcm^^ |54], although the properties 
of the stars that subsequently form under these conditions 
remain to be determined. Stacy & Bromm ll63ll have dis- 
cussed the impact of cosmic rays produced from Pop III 
stars, which would also be highly penetrating, although 
estimates of their production rate are uncertain. 

In gas that is no longer exposed to significant ionizing 
and dissociating flux and that has a chance to recombine, 
the relatively high residual electron fraction catalyzes 
molecule formation, particularly of HD, which can dra- 
matically enhance subsequent coolin g, pe rhaps reducing 
the characteristic star-formation mass i20ll21ll . Such a re- 
duction is in qualitative agreement with our models with 
lower values of K' , which is proportional to the temper- 
ature (eq. [T]i: such stars reach the main sequence (and 
associated higher ionizing luminosities) at lower proto- 
stellar masses. They also have smaller rates of accretion 
that they need to disrupt via disk photoevaporation. The 
precise mass range that is applicable for typical Pop III. 2 
stars requires more detailed numerical and analytic study, 
but initial estimates, based on the K' — 0.5 model, sug- 
gest that their masses are significantly smaller than those 
necessary to produce pair instability supernovae i60ll . 



THE TRANSITION TO POPULATION II 

Metal atoms and dust grains provide additional mecha- 
nisms for the gas to cool. For a given temperature, den- 
sity and ionization history there is some "critical metal- 
licity" at which these will start to influence star formation 
by dominating over H2 and HD cooling. 

The value of this critical metallicity for star formation 
is uncertain, with estimates ranging from r-^ lO^^Z© if 
the cooling is dominated by small d ust g rains that contain 
a significant fraction of the metals lUSll to ^ lO^^^Zg, if 
the cooling is dominated by the fine structure lines of C 
and O and there is negligible H2 |19]; if H2 cooling is 
included, Jappsen et al. ||64ll argue that there is no critical 
metallicity for gas-phase metals. 

A detailed understanding of the transition to Pop II 
star formation will require knowledge of the Pop III. 1 
and III. 2 IMFs, their mass loss during stellar evolution 
due to stellar winds and the metal and dust content of 
these winds, the metal and dust content of their super- 



nova ejecta, the mixing of wind and supernova ejected 
metals and dust into the surrounding (dense) gas that sub- 
sequently becomes gravitational unstable. It is fair to say 
that such a complicated problem, with so many uncertain 
inputs, remains far from solved. The general consensus 
is that, with additional coolants the gas can fragment into 
smaller pre-stellar cores that will tend to have smaller ac- 
cretion rates. Both of these effects should lead to an IMF 
peaked at smaller masses than that of Pop III stars. 



COMPARISON TO PRESENT-DAY STAR 
AND STAR CLUSTER FORMATION 

Initial Conditions 

The mean gas densities and velocity dispersions of 
Pop III.l minihalos and present-day star-forming gas 
clumps 116511 are quite similar The major difference is that 
in present-day gas, cooling by molecules and dust pro- 
ceeds down to ^ 10 K so that thermal pressure becomes 
negligible compared to other sources, such as magnetic 
and turbulent pressure. A magnetized, supersonically- 
turbulent, self-gravitating gas clump is prone to fragmen- 
tation: stochastic shock compressions produce a range of 
unstable core masses down to about the thermal Jeans 
mass. This process is likely to be responsible for shaping 
the present-day stellar IMF. Although Pop III. 1 cores are 
permeated by mild (sonic) turbulence, the compressions 
produced by these fluctuations are quite weak, and, cou- 
pled with the weak cooling properties of primordial com- 
position gas, means that fragmentation is suppressed. As 
a result, the IMF of Pop III. 1 stars is likely to be qual- 
itatively different from that of Pop II and Pop I stars: it 
should be set by the distribution of the entropy parameter 
K' and the rotation parameter /xep ■ Most numerical sim- 
ulations of Pop III. 1 star formation have not exhibited 
fragmentation in the minihalo down to their resolution 
limits, unless it occurs in a rotationally supported disk. 
Even this latter case appears to be unlikely to occur for 
realistic cosmological initial conditions. 

A second difference is that the sites of Pop III. 1 star 
formation are determined by dark matter potentials, i.e. 
the centers of minihalos. The cooling and settling of 
baryons mean that these typically dominate the mass 
densities on scales inside ~ 1 pc, but the dark matter still 
affects the gravitational potential of the larger scale core, 
which can affect the efficiency of feedback mechanisms. 
Dark matter annhilation in the centers of the halos could 
provide an important heating mechanism that could de- 
lay the onset of star formation Il33ll . 



Feedback Mechanisms 

One may ask how the feedback mechanisms of Pop 
III.l stars relate to those that operate in contempo- 
rary massive star formation. We note that the maximum 
mass attained in the fiducial model of McKee & Tan of 
Pop III.l star formation |15] is very similar to that in- 
ferred observationally in local massive star clusters (e.g. 
1,66.1 ). However, after decades of study, it remains unclear 
whether the maximum mass of stars forming today is set 
by feedback or instabilities in very massive stars [67]. 
We have argued that the maximum mass of primordial 
stars is set by feedback. The primary differences in the 
feedback processes then and now are: 

(1) Dust. In contemporary star- forming regions, dust 
destroys Lyman- a photons, eliminating them as a signif- 
icant pressure. On the other hand, the dust couples the 
pressure of the UV continuum radiation to the gas very 
effectively, and it remains to be determined whether this 
limits the final mass of the star; e.g., Yorke & Sonnhalter 
06811 find that it does, whereas Krumholz et al. |69] have 
not found evidence that it does. Dust also affects the evo- 
lution of H II regions, absorbing a significant fraction of 
the ionizing photons in dense H II regions, thereby reduc- 
ing the impact of H II region breakout and disk photoe- 
vaporation feedback. 

(2) Magnetic fields. In contemporary protostars, mag- 
netic fields drive powerful winds that drive away a sig- 
nificant fraction of the core out of which the star is form- 
ing 150(1. The cavities created by these winds allow radi- 
ation to escape from the vicinity of the protostar, signif- 
icantly reducing the radiation pressure |70]. As we have 
discussed, it is uncertain whether primordial protostars 
have such powerful magnetically-driven outflows |40]. 
However, even if they are present, we have argued that 
they would not be the dominant mechanism compared to 
disk photoevaporation in setting the final stellar mass. 

(3) Stellar temperatures and luminosities. Primordial 
stars were significantly hotter than contemporary stars, 
resulting in significantly greater ionizing luminosities. 
In addition, the accretion rates of primordial massive 
stars are much greater, at least initially, than those of 
contemporary massive stars li7 IJ . 
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